Classical Effects of Laser Pulse Duration on Strong-field Double Ionization 
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We use classical electron ensembles and the aligned-electron approximation to examine the effect 
of laser pulse duration on the dynamics of strong-field double ionization. We cover the range of 
intensities 10 14 — 10 16 W/cm 2 for the laser wavelength 780 nm. The classical scenario suggests that 
the highest rate of recollision occurs early in the pulse and promotes double ionization production 
in few-cycle pulses. In addition, the purely classical ensemble calculation predicts an exponentially 
decreasing recollision rate with each subsequent half cycle. We confirm the exponential behavior by 
trajectory back-analysis. 

PACS numbers: 32.80.Rm, 32.60.+i 



Detailed study of the classical picture of two-electron 
and multi-electron ejection from atoms and molecules is 
becoming important as short-pulse and high-field exper- 
imental work moves into the multi-electron regime 0. 
This is because it is well recognized that there is no re- 
alistic prospect of a quantum mechanical treatment of 
atomic and molecular multi-electron phenomena under 
the relevant experimental conditions, i.e., under highly 
non-perturbative fields with femtosecond and shorter 
time resolution. The exact solution of Schrodinger's 
time-dependent equation for helium double ionization in 
this regime by the Taylor group Q is accepted as the 
farthest that a comprehensive quantum treatment will 
go in the foreseeable future. At the same time, evidence 
has been accumulating that a purely classical scenario 
[jj 0] is consistent with many of the double ionization 
effects that are observed, and can be adapted further for 
multiple electrons. 

Simultaneous ejection of two electrons from atoms 
more complex than helium has been under experimental 
study for more than two decades [j| , and has been studied 
intensely in the decade since the discovery and confirma- 
tion [6| of a knee-like structure in the experimental ion 
count measurements in several inert gases, now accepted 
as implying anomalously high electron correlation. This 
signature has been reported for some molecules as well 
7J. Theoretical work based on a sequential ionization 
process involving independent electrons [j| is known to 
be inadequate to produce the knee signature. 

The prevailing picture of double ionization builds on 
the familiar electron rescattering scenario proposed for 
NSDI by Corkum . In this picture, one electron is ion- 
ized by the Keldysh mechanism |10| , tunneling out of the 
nucleus through a Coulomb barrier that has been lowered 
by the strong laser field. Then the laser field returns it in 
a purely classical fashion to the nucleus to collide with the 
bound electron and free it. A relatively direct classical 
analysis, without tunneling, of the effect can be based on 
solutions of Newtonian equations for individual electron 
trajectories obtained from large ensembles of 100,000 - 
1,000.000 members, which can be interpreted in a clear 
way |ll| . An example of a pair of energy trajectories 
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FIG. 1: (Color online) Plots of electron energy versus time 
measured in laser cycles, showing a pair of two-electron tra- 
jectories, with the trapezoidal laser field envelope fit) super- 
posed for reference. During laser turn-on the bound electrons 
exchange energy rapidly, and soon after laser turn-on is com- 
plete one (lighter red line) escapes and then returns in re- 
peated e-e recollisions until the inner electron (darker blue 
line) also has enough energy to escape the nucleus, and thus 
the pair become doubly ionized. Notice that after each rec- 
ollision, the energy of the inner electron is modified, and the 
outer electron carries different energies in each recollision. 



is shown in Fig. ^ where a sequence of rescatterings is 
clearly evident. Agreement between such classical and 
more elaborate partly or entirely quantum calculations 
EUll 13 E3 is remarkably good [IE El Q1 . 

Here we use the classical ensemble method with 
100,000 2-e trajectories and laser pulses with differ- 
ent numbers of laser cycles over a range of intensities, 
10 14 — 10 16 W/cm 2 . Each laser pulse has the experimen- 
tal wavelength 780 nm. Our laser field is taken to be a 
sine wave 19] with maximum field £ and a trapezoidal 
envelope function 1 > f(t) > with 2-cycle turn-on and 
turn-off, which is shown superposed on Fig. ^ For this 
report the number of cycles in the plateau ranges from 
to 10. The trapezoidal shape is chosen because it pro- 
vides a fixed plateau intensity and is known to suppress 
effects associated solely with turn-on pcj . 

To begin we show in Fig. [21 that in the domain of 
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FIG. 2: Variation of "knee" structure in laser pulses with 
different number of plateau cycles (iV). The probability of 
double ionization is the number of double ionization events 
divided by 100,000, the total number of 2-e trajectories used. 



interest here the knee signature of NSDI is preserved in 
a classical ensemble theory of the double ion count even 
for short pulses. We see that for intensities above the 
main knee region the double ionization yield is basically 
not dependent on the number of cycles in the pulse. We 
will see further evidence that in this regime the held is 
sufficiently strong to effect both electrons within a single 
cycle. 

As Fig. ^shows, while the pulse is turning on a rapidly 
repeated energy exchange occurs between the bound elec- 
trons, and this is what permits one to escape (with- 
out tunneling) over the Coulomb barrier, which becomes 
more and more suppressed as the field turns on. Our en- 
semble calculations 0) establish that when many cycles 
of laser pulse are available Fig. ^ is typical of the clas- 
sical process, i.e., the typical double-ionization event has 
been preceded by several or even many recollisions. This 
picture stands in two-way contrast to the "quantum" pic- 
ture in which tunneling is employed to release the first 
electron, and where double ionization is deemed to occur 
in the very first recollision cycle. 

There has been no way to choose between these two 
pictures of the double ionization timing. Experiments do 
not reveal how many cycles pass before the first electron 



is released. This uncertainty has be gun to be eliminated 
by the ability to use shorter pulses [21|, which can sup- 
press recollision effects requiring more than one cycle. 
These experiments show that double ionization does not 
require multiple recollisions, consistent with the "quan- 
tum" picture. The open question is whether the classi- 
cal theory, which has multiple-recollisions as typical, and 
doesn't invoke tunneling, agrees with this. The preser- 
vation of the knee signature in short pulses as shown in 
Fig. |21 confirms that NSDI can still proceed in the ab- 
sense of the typical multiple-recollisions. This positive 
answer prompts another question: how does this classi- 
cal picture incorporate both of these long and short pulse 
duration NSDI schemes? As we show here, the answer 
further confirms the attraction of a purely classical pic- 
ture. 

In order to accomplish classical double ionization the 
outer electron must find the inner electron in a suitable 
phase relation to its motion while going through the nu- 
clear region. We have already demonstrated [16j the effi- 
cacy of well-phased (what were called "slow-down" ) rec- 
ollisions. Of course, the outer electron may or may not 
return to the nucleus to make a recollision. If it returns 
to the nucleus, the recollision may take place in a condi- 
tion that is more or less favorable than in the previous 
encounter. Because of the energy exchange that initially 
helped free the outer electron, the inner electron is gener- 
ally to be found near the bottom of the nuclear potential 
well, where its bound oscillations are rapid and random. 
In a large ensemble it is reasonable that there will always 
be some in the right phase relation, even very early in the 
pulse. 

Our data confirms this qualitative picture, and pro- 
vides features not available in a quantum calculation, i.e., 
data based on specific recollision counts. This can pro- 
vide insights that are otherwise inaccessible. The first 
example is shown in Fig. [3J Here we show a new form of 
time record for pulses with 16 different intensities. We 
count the number of trajectories that fit the "recollid- 
ing" definition at 256 time instants, i.e., 32 times 
per cycle, for 8-cycle pulses at 16 different intensities. In 
these records the intensities are low enough that the en- 
tire recollision sub-ensemble is smoothly depleted during 
the double ionization over 8 cycles, and one can antici- 
pate a straightforward analysis, as follows. 

With each successful recollision some of the available 
population is removed. This means that the next recolli- 
sion has a fraction fewer candidates available for double 
ionization. A fractional decrease per half cycle is the 
same as predicting an exponential decrease in the num- 
ber of recollisions over time. To make a check of this pre- 
diction the classical approach suggests use of our "back- 
analysis" method 0] , but the number of successful dou- 
ble ionizations in any half cycle is very small even in a 
100,000-member ensemble, making statistical reliability 
difficult to achieve. However, the number of recollisions, 
successful or not, is much greater and can be used as a 
surrogate for the number of successful ones. Fig. 0]shows 
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FIG. 3: Number of recolliding electrons vs. time, counted at 8 times in each half cycle, for sixteen different intensities of 8-cycle 
pulses, from 4 x 10 13 to 10 15 W/cm 2 . The pulses with intensities in this range, which includes the main knee region, remove 
electron pairs smoothly and steadily, leading to the almost exactly exponential decrease of peak numbers through the duration 
of the pulse. 
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FIG. 4: (Color online) (a)-(d) Total number of recolliding 
electron trajectories in the overall ensemble counted at eight 
half cycles of an 8-cycle trapezoidal pulse with a 2-cycle on- 
ramp and 2-cycle off-ramp. Data is shown for four laser in- 
tensities in the NSDI knee region, and intermediate intensities 
produce results in the same range. The solid line in each plot 
corresponds to the best fit exponential function characterized 
with an exponential rate 7 for the data points. 



data extracted from Fig. Q|that confirms the exponential 
prediction and gives the exponential rates. 

One clear feature quickly seen in this new form of data 
is that the half cycle that typically has the most rec- 
ollisions, gateways to NSDI, is the very first one after 
the pulse reaches full strength. By itself this establishes 
that the classical picture does not neccessarily require 
many cycles of recollision in order to generate NSDI elec- 
trons. Since the more typical double ionization occurs af- 



ter several or more cycles, the absence of these multiplc- 
recollision gateways in short pulses reduces the produc- 
tion of double ionization as indicated in Fig. However, 
at only sightly higher intensities the picture changes. 

In Fig. OOthe record shows the same type of data shown 
in Fig. |21 but for a range of higher pulse intensities. 
The difference in effect is substantial, not only because 
the higher range of intensities promotes more effective 
double ionization, but because it also alters the nature 
of the process. The recollision begins to occur twice as 
rapidly. That is, it starts at the halfway point of the 
turn-on. Furthermore, most recollisions occur prior to 
peak intensity. The higher intensities in Fig. [S] merely 
emphasize this effect. One can see it beginning as early 
as in the third row of Fig. [3]where recollision starts prior 
to peak intensity. At high intensities, this effect coupled 
with the finding that the double ionization yield does not 
depend much on the pulse duration suggests that NSDI 
events can be completed within a single cycle. 

In summary, we have shown that NSDI production is 
predicted by the purely classical theory to occur at short 
times well in advance of the typical double ionization 
time. We see that the pulse begins to be strong enough 
to disturb a substantial fraction of the ensemble even be- 
fore the plateau regime is reached at the beginning of the 
second cycle. In fact, this distortion is due to the stronger 
effect of barrier suppression, which promotes one of the 
bound electrons to leave the nucleus, and allows recolli- 
sions to take place early in the pulse, which substantially 
enhances the NSDI production from a short pulse. 

This work was supported by DOE Grant DE-FG02- 
05ER15713. 
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FIG. 5: Number of recolliding electrons vs. time, counted at 8 times in each half cycle, for eight different intensities of 8-cycle 
pulses, from 2 x 10 15 to 9 x 10 15 W/cm 2 . The intensities in this range begin to be high enough to distort the NSDI process. 
The doubled peaks indicate that a major fraction of all members of the ensemble can be induced to leave the nuclear region in 
the same direction at the same time. 
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